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Abstract 

We make a systematic investigation on the two-body nonleptonic decays Be — s- D^^(^P, D^^!^V, 
by employing the perturbative QCD approach based on kx factorization, where P and V denote 
any hght pseudoscalar meson and vector meson, respectively. We predict the branching ratios 
and direct CP-asymmetries of these Be decays and also the transverse polarization fractions 
of Be — )• -^(^s)^ decays. It is found that the non-factorizable emission diagrams and annihi- 
lation type diagrams have remarkable effects on the physical observables in many channels, 
'■^ " especially the color-suppressed and annihilation-dominant decay modes. A possible large direct 

CP-violation is predicted in some channels; and a large transverse polarization contribution that 
can reach 50% ~ 70% is predicted in some of the Be — )■ D*s)^ decays. 
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I. INTRODUCTION 



The Be meson is the only quark-antiquark bound system (be) composed of both heavy 
quarks with different flavors, and are thus flavor asymmetric. It can decay only via weak 
interaction, since the two flavor asymmetric quarks (b and c) can not annihilate into 
gluons or photon via strong interaction or electromagnetic interaction. Because each of 
the two heavy quarks can decay individually, and they can also annihilate through weak 
interaction, meson has rich decay channels and provides a very good place to study 
nonleptonic weak decays of heavy mesons, to test the standard model and to search for 
any new physics signals 

Since the current running LHC collider will produce much more Be mesons than ever 
before, a lot of theoretical studies of the nonleptonic Be weak decays have been performed 
us,„g different approacl,es. For example. tl,e spectator-model h. the l.ghWrorrt quark 
model (LFQM) [3|, the relativistic constituent quark model (RCQM) [5|, the QCD 
factorization approach (QCDF) [6|, the Perturbative QCD approach (pQCD) 0jiol, and 
so on. Among the numerous decay channels, there is one category with only one charmed 
meson in the final states. They are rare decays, but with possible large direct CP asym- 
metry, since there are both penguin and tree diagrams involved. These decays have ever 
been studied in Ref. P] using the naive factorization approach. But they consider only 
the contribution of current-current operators at the tree level, and thus no direct CP 
asymmetry is predicted. They also have difficulty to predict those pure penguin type 
or annihilation dominant type decays, such as Be — > D^(f), D^K^, D'^(j). Ref. js] dis- 
cussed some semileptonic and nonleptonic Be weak decays and CP-violating asymmetries 
by using RCQM model based on the Bethe-Salpeter formalism. They do not include 
the contributions of annihilation type diagrams, either. Since the annihilation type con- 



tributions are found to be important in the B meson non-leptonic decays Uj and also 



significant in the Be decays [l2j, one needs further study these channels carefully. 

(*) 

In this paper, we calculate all the processes of a Be meson decays to one D^^^ meson and 
one light pseudoscalar meson (P) or vector meson (V) in pQCD approach. It is well-known 
that Be meson is a nonrelativistic heavy quarkonium system. Thus the two quarks in the 
Be meson are both at rest and non-relativistic. Since the charm quark in the final state 
D meson is almost at coUinear state, a hard gluon is needed to transfer large momentum 
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to the spectator charm quark. In the leading order of rric/mB,. ~ 0.2 expansion, the 
: actorization theorem is apphcable to the Be system similar to the situation of B meson 
isj. Utilizing the kr factorization instead of collinear factorization, this approach is free 
of endpoint singularity. Thus the diagrams including factorizable, nonfactorizable and 
annihilation type, are all calculable. It has been tested in the study of charmless B meson 
decays successfully [l^, especially for the direct CP asymmetries 15 1. For the charmed 



decays of B meson, it is also demonstrated to be applicable in the leading order of the 
mo/f^B expansion jl6[j2l|. 

Our paper is organized as follows: We review the pQCD factorization approach and 
then perform the perturbative calculations for these considered decay channels in SecJITl 
The numerical results and discussions on the observables are given in Sec JIIII The fi- 
nal section is devoted to our conclusions. Some details related functions and the decay 
amplitudes are given in the Appendix. 

II. THEORETICAL FRAMEWORK 

For the charmed B^. decays we considered, the weak effective Hamiltonian T-ieff for 
h — q'{q' = d, s) transition can be written as [22 1 

^ 10 

* q=u,c i=3 

with the Cabbibo-Kobayashi-Maskawa (CKM) matrix element = VqqiV*^^. Oi{fi) and 
Ci{fi) are the effective four quark operators and their QCD corrected Wilson coefficients. 



respectively. Their expressions can be found easily for example in Ref. [22 1. 

With these quark level weak operators, the hardest work is left for the matrix ele- 
ment calculation between hadronic states {DM\Heff\Bc) ■ Since both perturbative and 
non-perturbative QCD are involved, the factorization theorem is required to make the 



calculation meaningful. The perturbative QCD approach 1^ is one of the methods to 
deal with hadronic B decays based on factorization. At zero recoil of D meson in the 
semi-leptonic Be decay, both c and b quark can be described by heavy quark effective 
theory. However, when the D meson is at maximum recoil, which is the case of two body 
non-leptonic B^ decay, the final state mesons at the rest frame of B^ meson are collinear, 
so as to the constituent quarks (c and other light quarks) inside. Since the spectator 
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c quark in the Be meson is almost at rest, a hard gluon is then needed to transform it 
into a colhnear object in the final state meson. This makes the perturbative calculations 
into a 6-quark interaction. In this colhnear factorization calculation, endpoint singular- 
ity usually appears in some of the diagrams. The QCD factorization approach ^] just 
parameterize those diagrams with singularity as free parameters; while in the so-called 
soft-collinear effective theory 2J], people separate these incalculable part to an unknown 
matrix element. In our pQCD approach, we studied these singularity and found that they 
arise from the endpoint where longitudinal momentum is small. Therefore, the transverse 
momentum of quarks is no longer negligible. If one pick back the transverse momentum, 
the result is finite. 

Because the intrinsic transverse momentum of quarks is smaller than the b quark mass 
scale, therefore we have one more scale than the usual colhnear factorization. Additional 
double logarithms appear at the perturbative QCD calculations. These large logarithms 
will spoil the perturbation expansion, thus a resummation is required. This has been 
done to give the so called Sudakov form factors [25] . The single logarithm between the W 
boson mass scale and the factorization scale t in pQCD approach has been absorbed into 
the Wilson coefficients of four quark operators. The decay amplitude is then factorized 
into the convolution of the hard subamplitude, the Wilson coefficient and the Sudakov 
factor with the meson wave functions, all of which are well-defined and gauge invariant. 
Therefore, the three-scale factorization formula for exclusive nonleptonic B meson decays 
is then written as 



C{t) (g) H{x, t) (g) ® exp 



.(P,6)-2 f ^7,(«s(/x)) 
Ji/b 



(2) 



h/b A* 

where C{t) are the corresponding Wilson coefficients. The Sudakov evolution 
exp[— s(P, 6)] I25I are from the resummation of double logarithms ln^(P6), with P denot- 
ing the dominant light-cone component of meson momentum. 7^ = —Us/t^ is the quark 
anomalous dimension in axial gauge. All non-perturbative components are organized in 
the form of hadron wave functions which can be extracted from experimental data or 

other non-perturbative method. Since non-perturbative dynamics has been factored out, 
one can evaluate all possible Feynman diagrams for the six-quark amplitude straight- 
forwardly, which include both traditional factorizable and so-called "non-factorizable" 
contributions. Factorizable and non-factorizable annihilation type diagrams are also cal- 
culable without endpoint singularity. 



The meson wave function, which describes hadronization of the quark and anti-quark 
inside the meson, is independent of the specific processes. Using the wave functions 
determined from other well measured processes, one can make quantitative predictions 



here. For the hght pseudoscalar meson, its wave function can be defined as 
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$(P,x,0 = ^l,[r^i{x)+mo<Pl{x)+imo{iiij - l)0^(x)], (3) 

where P is the momentum of the light meson, and x is the momentum fraction of the 
quark (or anti-quark) inside the meson. When the momentum fraction of the quark (anti- 
quark) is set to be x, the parameter ^ should be chosen as +1(— 1). The distribution 
amplitudes 0p(x),0p(a;) and 0p(x) are given in Appendix [Cl 

For the light vector mesons, both longitudes (L) and transverse (T) polarizations are 
involved. Their wave functions are written as [3] 

$^(X) = -^{Myfy'^Mx) + ivP<Pv{^) + Mv<PU^)}ap. 
V ^jVc 

$^(x) = -^{M^/f 0^(x) + f^pfvi^) + iMve^,p,i,e*T^nPv''ry{x)}o.p, (4) 

where e^*-"^^ denotes the longitudinal (transverse) polarization vector. And convention 
g0i23 = X is adopted for the Levi-Civita tensor. The distributions amplitudes are also 
presented in Appendix [Cl 

Consisting of two heavy quarks (b,c), the Be meson is usually treated as a heavy 

aiarkonium system. In the non-relativistic limit, the Be wave function can be written as 
1 

^bM) = ^[(Z' + Mb.)i,5{x - rj], (5) 

with Te = rrie/MB^- Here, we only consider one of the dominant Lorentz structure, and 
neglect another contribution in our calculation j27| . 

In the heavy quark limit, the two-particle light-cone distribution amplitudes of 
D(s)/-D('g) meson are defined as [21| 

(D(,)(P2)|g.(^)c/3(0)|0) = /'dxe"^-^[75(yp2 + mB(j0B(^,(a;,6)],;5, 

{DU{P2)\q.{z)cp{m) = -^^'rfa;e"^^-1/(A + ^i>tJ%)^^'^)]"/^- 
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We use the following relations derived from HQET to determine fo* 28 1 




(7) 



For the D\*! meson wave function, we adopt the same model as of the B meson [8| 



with shape parameters ud = 0.6 for D/D* meson and ojds = 0.8 for Dg/D* meson. Here, 
a larger ud^ parameter than ud characterize the fact that the s quark in Dg meson carries 
a larger momentum fraction than the light quark (u,d) in the D meson. 

At leading order, there are eight types of diagrams that may contribute to the 
Be — > i5|*jp, -D(*jv^ decays as illustrated in Fig{Tl The first line are the emission type 
diagrams, with the first two contributing to the usual form factor; the last two so-called 
non-factorizable diagrams. The second line are the annihilation type diagrams, with the 
first two factorizable; the last two non-factorizable. 



(x, b) = Nd^^^ [x{l - x)]^ exp 
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A. Amplitudes for Be — t- ^(s)^ decays 

We mark LL, LR, SP to denote the contributions from {V-A){V -A), {V-A){V + A) 
and {S — P){S + P) operates, respectively. The amphtudes from factorizable diagrams 
(a) and (b) in FigJT] following: 

^^^^ = 2^CffBfpirM^ dx^j^ b,b2dbidb2Mx2,b2) X 

{[(1 - 2rD)x2 + {td - 2)n]as{ta)he{ae, /3a, bi, b2)St{x2) exp[-S'afe(ta)] 
-{td - 2)rD(a;i - l)as{tb)he{ae, /3b, &2, bi)St{xi) exp[-Sab{tb)], (9) 

where td = itid/Mb, h = nib/MB', Cp = 4/3 is a color factor; fp is the decay constant 
of pseudoscalar meson (P). The factorization scales ta^b are chosen as the maximal vir- 
tuality of internal particles in the hard amplitude, in order to suppress the higher order 
corrections {29 1. The function he are displayed in the Appendix iBl The factor St{x) is the 
jet function from the threshold resummation, whose definitions can be found in [8J. The 
terms proportional to r|, have been neglected for small values. We can calculate the form 
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factor from eq.Q if we take away the Wilson coefficients and fp. For the {V — A){V + A) 
operates, we have J^^^ = —J^^^ since only axial-vector current contributes to the pseu- 
doscalar meson production. For the [S — P)(5' + P) operates the formula is different: 

Jf^ = -A^CffsfpTTM^s dx^j^ hb2dhdb2M^2,b2) X 

{[rD{4:rb-X2 - 1) - n + 2]as{ta)he{ae, (3a,bi,b2)St{x2) exp[-Sab(ta)] 

+ [r£i(2 - Axi) + xi]as(tb)he{ae, /3b, &2, bi)St{xi) exp[-Sab{tb)]. (10) 

For the nonfactorizable emission diagram (c) and (d), the decay amplitudes of three 
types operates are: 

8 

= -CpfB-^M^ dx2dxs b2b3db2db3(pD{x2,b2)(pi{x3) X 
<J Jo Jo 

{[rz)(l - xi - X2) + X1 + X3- l]asiQK{(3c, a^, 63, exp[-Scd{tc)] - 
[rz)(l -xi- X2) + 2x1 + xa - xs - l\as{td)K{Pd, ol^-, h, h) exp[-5'cd(td)]}, 



>i ^00 



-^e = -CpfBTTM^rpil + rD) I dx2dx3 I b2b^db2db:i(pD{x2M) ^ 



JO 

Pi 



{[{Xi + X3 - 1 + rz5(2xi + X2 + X3 - 2))0p(x3) + 

(xi + X3 - 1 + r£,(x3 - X2))(tPp{x'i)]as{tc)K{l3c, "e, &3, h) exp[-S'cd(tc)] 

-[(Xi - X3 + rz)(2Xi + X2 - X3 - l))0p(x3) + (X3 - Xi + 

roix^ + X2- I))(f)p{x3)]as{td)he{f3d, Oe, &3, h) exp[-Scd{td)]}, (12) 



8 C 

Mf^ = -CpfBirM^ dx2dx3 b2b3db2db3(t)D{x2M)(t>j'{xz) x 
<J Jo Jo 

{[rnixi + X2 - 1) - 2xi - X2 - X3 + 2] 
asitc)he{/3c, Oe, h, h) exp[-Scditc)] - 

[xs - xi - rr){l - xi - X2)]as{td)K{(3d, a^, &3, ^2) exp[-Scditd)]}, (13) 

where rp = rriQ /Mp, with as the chiral mass of the pseudoscalar meson P. 

For the factorizable emission diagram (e) and (f), we keep the mass of the c-quark in 
D meson, while the mass of the light quark is neglected. The amplitudes are given as 
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follows: 

/•l roc 

Jo Jo 

{[(pj,{x3){x3 - 2rDrc) + rp[(f>^{x3){2rD{x3 + 1) - re) + (pUxs) 

{Vc + 2rD{x3 - l))]]as{te)he{aa, Pe, &2, ^s) exp[- Sef{te)]St{x3) - 

[x2<Pi{x3) + 2rprD{x2 + 1)0^(0:3)] 

as{tf)he{aa, 13 f, 63, 62) exp[-S'e/(i/)]S't(x2)}, (14) 
Jo Jo 

{[-(pp{x3){2rr) - Tc) + rp[(pp{x3){4rcrD - X3) + (j)p{x3)x3W 

as{te)he{aa, /3e,b2,b3) exp[-Sef{te)]St{x3) - [x2rD4>p{x3) + 

2rp(l)p(x3)]as(tf)he{aa, Pf, 63, ^2) exp[-Sef(tf)]St(x2)}; (15) 
and that of the nonfactorizable annihilation diagram (g) and (h) are 

M^^ = -^CpfBT^Mp j dx2dx3 / bib2dbidb2(f)D{x2,b2) X 
-J Jo Jo 

{[4>iix3){rc -xi + X2) + rprD[4'pix3){x2 - X3) + 

(f)p{x3){irc - 2xi +X2 + X3)]]as{tg)he{/3g, 61, 62) e^[-Sgh{tg)] 

+ [-(piix3){rb + Xi+ X3 - 1) + rprD[ix2 - X3)(f)'§,{x3) - 0p(x3) 

(4r6 + 2x1 ^X2^X3- 2y^a s{th)he{Ph, oca-, bi, 62) exp[-Sgh{th)]}, (16) 

= -CfJbT^M^ / cix2CiX3 / bib2dbidb2(pD{x2,b2) X 

•J Jo Jo 

{[-(t)pix3)rD{rc + xi- X2) + rp[-(f)p{x3){-rc - xi + X3) 
+<Ppix3){rc + xi- X3)]]asitg)heiPg, aa, bi, 62) exp[-Sghitg)] + 
[-(pp{x3)rD{-rb + Xi + X2-1)+ rp[{-rb + X1 + X3-I) 
{(ppixs) + (j)p{x3))]]as{th)he{Ph, aa, bi, 62) exp[-Sgh{th)]}, (17) 

8 /"^ f°° 

J^a^ = --CpfsT^Mp / dX2dX3 / &i&2(^fei(^&20D(a:^2, ^2) X 

<J Jo Jo 

{[-0p(x3)(a;i - X3 - Tc) + rprD[-(f)p{x3){x2 - X3) 

+(f)p{x3){irc - 2xi + X2 + X3)]]as{tg)he{/3g, aa, bi, 62) e^[-Sgh{tg)] + 

[-(j)j,{x3){rb + xi +X2-1)+ rprz)[(-4r5 - 2xi - X2 - X3 + 2)0^(^3) 
-{x2 - X3)(pp{x3))]]as{th)he{Ph, Ola, bi, 62) exp[- S gh{t h)]} , (18) 
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With the functions obtained in the above, the total decay amphtudes of 10 decay channels 
for the Be — > D(^s)P can be given by 

+(^6 + + + + -FD], (19) 



LL 



+ (-C5 + \c7)M^e'' + (-C4 - - ClO - ^Cg) + 

5 1 3 1 LL 

(ClO + gCg - -Cs - C4 - -C7 - -Cs)J^e 

+(-C6 - ^Cs + \Cs + ^^7) - (C5 + C7)A^^^ 

+(-^6 - \c, -Cs- lc,)J^r], (20) 

-6[(2C4 + Cs + ^Cio - ic9)Mf + (Cs + C9)A^f 
+(C5 - ^C7)Mf ^ + (C5 + C7)A4^^ + (C4 + ^Cs + Cio + 
^Cg) J-f^ + (^Ca + ^C4 + i(C9 - Cio))-F,^^ + (2C5 + ^Ce 
+1^7 + iCs) 7-,^^ + (Ce + ^C5 - ^Cs - ^C7) JT 
+(C6 + ^C5 + C8 + ^C7)^fn, (21) 

A{B, ^ L'+Ty,) = -6[(C4 - ^Cio)A^f + (Ce - ^Cs)^^ + (C3 + ^C4 

-1^9 - iCio)^^^ + (C5 + ^Ce - ic7 - iCs) J-i'^], (22) 

v^^(S, ^ D+ry,) = e:[«2-^i^^ + C2A^f^] - ea(2C4 + ^Cio)Me^^ + 

(2(^6 + ^C8)Alf + (2^3 + ^C4 + ^Cg + ^Cio)^e" 

+(2C5 + ^Ce + ic7 + lCs)J'^% (23) 
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+(C4 + + Cio + ^Cg) J-f ^ + (Ce - ^C8)7Wf^ + 

^(2(^3 + C4) - (Cg + Cio)) J-f ^ + (C5 + - \Cj - ^Q) J-f ^ 



(24) 

+(C3 + C9)M^^ + (C5 - ^CO-Me^^ + (C5 + C7)-Mf^ 

+ (C4 + ^C^3 + Cio + ^C^9)-^a''^ + 

-^C^7)-^f'' + (Ce + ^C^5 + + ^Cr) J-^], (25) 

^(i?,^D>°) = C[a2^e^^ + C2-M^^]-ei[(^(3C79 + Cio)-F,^^ 

+i(3C7 + Cs)) J-^^ + \c,oM^,^ + ^CsA^f ], (26) 

with the CKM matrix element = VidV*^ and ^[ = VisV^l{i = u,c,t). The combinations 
Wilson coefficients ai = C2 + Ci/3 and 02 = Ci + C2/S. The total decay amplitude of 
A{Bc ^ D^K+) and ^(5^ ^ D+K°) can be obtained from ([19]) and ([25D, respectively, 
with the following replacement: 

A{B,^D'K+) = A{B, ^ D%+)l^K,i^^^i, 

A{B, ^ D+K') = A{B, ^ DtK'>)\n,^D,i,^^,. (27) 

It should be noticed that, in ( 12T|) . ( I22l) . ( |23|) and (El]), the decay amplitudes are for the 
mixing basis of {riq,r]s). For the physical state {rj,!]'), the decay amplitudes are 

AiB^ D+ri) = AiB^ D+r]q) cos(p- A{B^ D+r]^) sincf), 

A{Bc D+T]') = AiBc D^rig) sin0 + A{B^ D+r]s) coscp, 

A{B, ^ D+7]) = A{B, ^ D+7]g) COS0 - AiB, ^ D+r/,) sin0, 

AiB, Dtv) = A{B^ D+7]g) sin + A{B, D+r/,) cos 0, (28) 
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where = 39.3° is the mixing angle between the two states. 

f T] \ _ f cos(j) -sm(f)\ f r]q\ 

\r]'j~\sm(j) COS0 J \Vs J ' ^ ' 

B. Amplitudes for Be — )■ ^(s)^ decays 

In Be —7- -D(s)l/ decays, the vector meson is longitudinally polarized. In the leading 
power contribution, the formula of each Feynman diagram for — ?► -0(^)1^ is similar to 
that of the Be — )■ -D(s)P modes, but with the replacements 

fp fv, rp rv, (j)v, 0P ^ -0y> 0p ^ (Pv- (30) 

The total decay amplitude for Be — > -0(^)1^ can be obtained through the substitutions in 

TT ^ p, K ^ K*, r]g^ w, 7]s 0. (31) 

C. Amplitudes for Be — ^ ^(s)^ decays 

The decay amplitude of Be — )■ D*^-^V can be decomposed into 

A{eD* ,ev)=A^ + A^eD*T ■ evT + M^e^/jp^n^w'^e^.^e^^ (32) 

where en'ri^VT) is the transverse polarization vector for D*(y) meson. corresponds 
to the contributions of longitudinal polarization; A^ and A'^ corresponds to the contri- 
butions of normal and transverse polarization, respectively. The factorization formulae 
for the longitudinal, normal and transverse polarizations are all listed in Appendix |X1 
There are also 10 channels for Be — t- D'^s)^ decay modes. We can obtain the total decay 
amplitudes from those in Be — )■ -D(s)^ with replacing by D'^g)- 

D. Amplitudes for Be — ^ ^\s)^ decays 

For Be — > Z^^g^P, only the longitudinal polarization of D*^-^ will contribute. We can 
obtain their amplitudes from the longitudinal polarization amplitudes for the Be — )■ D*^^s)^ 
decays with the following replacement in the distribution amplitudes: 

fv fp, rv rp, (f)v 0p, 0y ^ 01- ^ 0p- (33) 
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TABLE I: The decay constants and the hadronic meson wave function parameters taken from 
the hght-cone sum rules 33 1. 



The decay constants (MeV) 



489 



fD 



206.7 ±8.9 



257.5 ±6.1 



131 



fK 



160 



fp 



209 



165 



f. 



195 



T 



145 



f(p 



231 



J 6 



200 



Values of Gegenbauer moments 





vr 


K 


Vq 


Vs 


af 




0.17 








0.25 


0.115 


0.115 


0.115 




-0.015 


-0.015 


-0.015 


-0.015 




P 


LO 




K* 


4- 








0.03 


^2 


0.15 


0.15 


0.18 


0.11 










0.04 




0.14 


0.14 


0.14 


0.10 



fK- 



217 



Jk* 



185 



In fact, the — t- D*^s)^ decays amplitude are the same as the B^. D{s)P ones only at 
leading power under the hierachy Mb^ ^ ^ ^qcd- An explicit derivation shows 

that the difference between the two kinds of channels occurs at (9(r^(,)) and at the twist-3 
level in eq.fl9|)-eq. (fT8|) . 



III. NUMERICAL RESULTS AND DISCUSSIONS 



The numerical results of our calculations depend on a set of input parameters. We 
list the decay constants of various mesons and parameters of hadronic wave functions in 
Table [H For 77 — 77' system, the decay constants fg a nd fs in the quark-flavor basis have 



been extracted from various related experiments 30|, |31| 



/, = (1.07 ± 0.02)A, /, = (1.34 ± 0.06)/^. 



(34) 



For the CKM matrix elements, the quark masses etc., we adopt the results from 32 1 



\Vub\ 

\Vus\ 

rrir. 



(3.89 ± 0.44) X 10-^ \Vud\ = 0.97425, \Vcb\ = 0.0406, = 0.23 

0.2252, \Vcs\ = 1.023, 7 = (73^ 



-22 \o 

-25J ' 



1.27GeV, rrib = 4.2GeV, 



m. 



K 



1.6GeV, 7n[|' 



1.07GeV, 



771" = 1.4GeV, 



Tng^ = 1.92GeV, A^^cd = 0.112GeV. (35) 
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For the considered Be Di^g-^P^ Be — )■ D*^^P and Be -^(s)^ decays, the branching 
ratios BTZ and the direct CP asymmetry A'^p for a given mode can be written as 

'^'^~32nMB^ ^^^^ 
where the decay amphtudes A have been given exphcitly in Sec. [Ill for each channel 
A is the corresponding charge conjugate decay amphtude, which can be obtained by 
conjugating the CKM matrix elements in A. 

Our numerical results of CP averaged branching ratios and direct CP asymmetries for 
Be — )■ -D(s)P and Be — t- -D(s)^ decays are listed in Tables Hi] and IIIIl respectively. The 
dominant topologies contributing to these decays are also indicated through the symbols 
T (color-allowed tree), C (color-suppressed tree), P (penguin) and A (annihilation). The 
first theoretical error in all our tables is referred to the Z)(5) meson wave function: (1) The 
shape parameter Uo = 0.60±0.05 for D / D* meson and u)Ds = 0.80±0.05 for Dg/D* meson; 
(2) The decay constant fo = (206.7 ± 8.9)MeV for D meson and fo, = (257.5 ± 6.1)MeV 
for Dg meson. The second error is from the combined uncertainty in the CKM matrix 
elements Vub and the angle of unitarity triangle 7, which are given in eq. (l35|) . The third 
error arises from the hard scale t varying from 0.75t to 1.25t, which characterizing the size 
of next-to-leading order QCD contributions. Most of the branching ratios are sensitive to 
the hadronic parameters and the CKM matrix elements. The CP asymmetry parameter is 
only sensitive to the next-to-leading order contributions, since the uncertainty of hadronic 
parameters are mostly canceled by the ratios. 

We also cite theoretical results for the relevant decays evaluated in LFQM model 
and RCQM model jsl to make a comparison in Tables [TTl and IIIIl Our pQCD results are 
generally close to RCQM results but differ substantially from the ones obtained by LFQM. 
This is due to the fact that LFQM used a smaller form factors F^^^^^q^ = 0) = 0.086 
at maximum recoil, which is rather smaller than other model predictions [4] and also 



another covariant LFQM results 3J]. In fact, these model calculations all give consistent 
form factors at the zero coil region, considering only soft contributions by the overlap 
between the initial and final state meson wave functions, which is good at the zero recoil 
region. At the maximum-recoil region, which is the case for non-leptonic B decays, the 
soft contribution is suppressed, since a hard gluon is needed, as discussed in the previous 
section. Furthermore LFQM only consider the contribution of current-current operators at 
the tree level, therefore they cannot give predictions for those modes without tree diagram 
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TABLE II: CP averaged branching ratios and direct CP asymmetries for B(. — >■ -D(s)P decays, 
together with results from RCQM and LFQM. 



channels 



Be 

Be 
Be 

Be 
Be 
Be 
Be 
Be 
Be 
Be 



D + TT^ 

D+K^ 

D + T] 
D + TT^ 

D+ri' 



Class 



T 
C,A 
A,P 
A,P 
C,A 
C,A 
C,P 
A,P 
A,P 
A,P 



This work 



+3.1+6.0+0.8 
3.5-5.6-0.6 
+0.24+0.55+0.06 
-0.16-0.41-0.01 
+17.2+2.2+5.4 
-9.1-1.7-3.6 
15.6+0.3+7.4 
12.3-0.3-4.6 
n Q9+0. 15+0.21+0.03 
'-'•^^-0.15-0.25-0.00 
0.12+0.16+0.06 
0.10-0.20-0.03 



26.7] 
0.82^ 

46.r 



o.gr 



0.41 



+0.04+0.01+0.02 



0.04-0.02-0.02 
9 1+0.9+0.3+0.3 
^•-•--0.6-0.3-0.2 
1.7+0.5+3.3 
1.8-0.6-1.2 



17.3^ 



51.0 



+4.9+0.4- 
-5.4-0.3- 



-6.7 
-3.5 



RCQM^ 



22.9 
2.1 
44.5 
49.3 



1.9 



LFQM 



4.3 
0.067 
0.35 

0.087 
0.048 
0.0067 

0.009 
0.0048 



This work 



-41.2 

2.3 
-34.8 

2.3 



+4.5^ 



11.1+0.8 
■4.6-7.8-1.2 
1+6.3+1.4+15.0 
3.0-0.8-18.8 
-4.9+7.4+1.8 
2.6-3.7-1.3 
,+0.4+0.9+0.0 
0.2-0.5-0.0 
An O+0.0+18.4+15.6 
^U.O_2 g_i4 0-13.5 

_-,4 n+0-6+4.6+15.9 
-•-^•^-1.5-5.2-11.9 
/Ifi 7+1.4+6.3+2.5 

' -1.4-11.8-2.8 
c-A Q+6. 9+5.3+0.0 
0^.0_7 2-8.0-0.3 
9 o+O.O+0. 4+1.1 
^•"-0.1-0.7-1.2 
-1 -1+0.1+0.2+0.7 
-•-•^-0.0-0.2-0.6 



RCQM 



6.5 
-1.9 
-4.6 
-0.8 



13.3 



"■we use the results of decay -widths in [Sj , but we take tb^ — 0.453ps to estimate the branching ratio 



TABLE III: CP averaged branching ratios and direct CP asymmetries for Be — )• D(^g)V decays, 
together with results from RCQM and LFQM. 



BTZjW 



-7\ 



channels 


Class 


This work 


RCQM 


LFQM 


This work 


RCQM 


Be- 


+ D0 




T 


fifi 9+7.6+16.0+1.6 
""•■^-7.6-14.1-1.3 


60.0 


13 


9/1 C-+2.6+5. 3+0.3 
^^•^-0.4-3.2-0.8 


3.8 


Be- 




P° 


C,A 


-1 4+0.1+0.5+0.1 
-"-•^-0.2-0.5-0.2 


3.9 


0.2 


7Q O+0.3+11.2+3.4 
' ^•"-5.8-19.6-10.7 


-3.0 


Be- 


+ Z)0 


K*+ 


A,P 


9C- Q+2.7+0.9+1.5 
^-^•^-3.0-0.8-0.8 


34.7 


0.68 


an 9+1.8+15.1+0.7 
""•^-0.6-6.5-0.0 


-6.2 


Be- 


+ Z)^ 




A,P 


-IQ -1+3.3+0.1+0.7 
-•-^•-•--2.5-0.0-0.7 


28.8 




Q C-+0. 0+0.5+0. 5 
"^•'-'-0. 1-0.8-0.3 


-0.8 


Be- 


^D^ 




C,A 


-I Q+0. 3+0.5+0.0 
-•-•^-0.3-0.6-0.0 




0.15 


q /.+3.9+1.3+13.4 
"^•"-1. 2-1. 6-10.7 




Be- 


+ Z)^ 




P 


"•"'-'"-0.001-0.0-0.001 










Be- 


+ Z)^ 


P° 


C,P 


n Qc+O.lO+O.02+0.04 
"•^^-0.09-0.01-0.04 






c-n 9+1.0+5.9+2.5 
"-•"•^-1.1-11.9-3.2 




Be- 


+ Z)^ 


-^*o 


A,P 


-1 /i+0. 2+0.0+0.1 
-•-•^-0.2-0.1-0.1 


1.0 




a^ 5+4.5 
"-•-•"-0.3-14.2-3.6 


13.3 


Be- 


^D^ 




CP 


rv Q-i +0.03+0.07+0.07 
"•"Jj--0.03-0.07-0.05 




0.016 


AA Q+0.8+17.1+10.3 
^^•^-1. 6-14.9-13. 6 




Be- 


+ Z)^ 




A,P 


97 n+4.8+0.1+2.0 
^' •"-1.2-0.0-0.4 


15.7 


0.0048 


q q+0. + 0.4+0. 3 
"^•"^-0. 3-0.8-0. 4 


-0.8 



contributions like Be — ?■ D^K^ and Be — ?■ D^K^. For the color suppressed decays (C), our 
predictions differ from the ones of RCQM, since in these modes, the contributions from 
the non-factorizable emission diagram and annihilation diagram dominated the branching 
ratio, which are not calculable in RCQM. 

Our numerical results of the CP averaged branching ratios and direct CP asymmetries 
for Be — 7- D*^s)-^ decays are listed in Table HVl together with the RCQM model predictions. 
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Again, our results are similar with RCQM model for the tree dominant mode (T). But 
for the annihilation dominant and penguin dominant modes (A,P), the branching ratios 
obtained in the RCQM are one order of magnitude smaller than ours. The reason is that 
these decay amplitudes are governed by the QCD penguin parameters 04 and ag in the 
combination 04 + Ra^ 35| in the factorization hypothesis. The coefficient R arises from 
the penguin operator Oq, where R > for B PP, R = for PV and VV final states, 
and -R < for B — )■ VP, the second meson in the final states is the one emitted from 
vacuum. Therefore, the branching ratios of various type decays have the following pattern 
in the factorization approach 

Bn{B, DP) > Bn{B, DV) ~ BTZ{B^ D*V) > BTZ{B^ D*P) (37) 

as a consequence of the interference between the 04 and ag penguin terms. In the contrary, 
we have additional non-factorization contributions and large annihilation type contribu- 
tions in the pQCD approach, which spoils the relation in eq. flS7|) . 

As expected, the annihilation type diagrams give large contributions in the B^. meson 
decays, because the annihilation type diagram contributions are enhanced by the CKM 
factor V*fycq 0, 3- For the 6 -)■ (i process, | | = 2.5; For the 6 s process, 

V* V 

= 47. The annihilation diagram contributions are the dominant contribution in 
some h ^ s processes. Therefore, we have the ratio relation ^^^^^(*)+l^(*)o^ ~ 1 for 
these two annihilation dominant b —> s transition processes. 

For the color suppressed decays (C), our predictions differ from the ones of RCQM, 
since in these modes, the contributions from the non-factorizable emission diagram and an- 
nihilation diagram dominated the branching ratio, which are not calculable in RCQM. For 
example in decays B^ — ?► D'^*^^(7i^,r],r]', p^,uj), the non-factorizable contribution, which 
is proportional to the large Wilson coefficient C2, is the dominant contribution. In fact, 
the annihilation diagrams can also give relatively large contributions for the enhance- 
ment by CKM factor. We also find that the twist-3 distribution amplitudes play an 
important role, especially in the factorizable annihilation diagrams. As stated in sec- 
tion III Dt the Be — > DP{V) decay amplitudes are different from Be — )■ D*P{V) ones 
only at twist-3 level. The numerical results show that the contributions from factorizable 
annihilation diagrams have an opposite sign between the two type channels. For exam- 
ple, this results in a constructive interference between non-factorizable emission diagrams 
and factorizable annihilation diagrams for B^ D*^tt^, but a destructive interference 
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TABLE IV: CP averaged branching ratios and direct CP asymmetries for Be ^(s)-^ decays, 
together with results from RCQM. 





B7^(lo-7) 




channels 


Class 


This work 


RCQM 


This work 


RCQM 


Be D*^Tr+ 

Be D^tt" 
Be D*^K+ 
Be D*+K^ 
Be -+ D*+ri 
Be D*+ri' 
Be L»*+7r° 
Be D*+K^ 
Be ^ D*+r^ 
Be ^ Dl+r]' 


T 

C,A 
A,P 
A,P 
C,A 
C,A 
CP 

A,P 
A,P 
A,P 


10 0+2.0+4.1+0.4 

10. 0_2 0-3.5-0.5 
1 q+0.4+0.2+0.0 
^•■^-0.3-0. 3-0.0 
7q c;+31.0+0. 8+0.7 
(0.0_23.4-l. 1-0.4 

77 0+25.4+0.2+7.2 
' ' •°-24.0-0.2-5.2 
n q/i+0. 14+0. 19+0. 04 
^•"-•^-0.09-0.15-0.00 
n 1 tr+0. 08+0. 08+0. 03 
"-'•■'-^-0.05-0.06-0.01 
n. O7+0.02+0.03+0.01 
"-'•^'-0.04-0.02-0.02 
1 /.+0.2+0. 1+0.2 
^•"-0.3-0.1-0.1 
7+5.3+0.3+0.1 
-•^"J^' -4.0-0.2-0.3 
14 4+6.6+0.1+0.5 
-'^^•^-4.6-0.1-0.6 


19.6 
0.66 
4.9 
2.8 

0.21 


f>4 n+12.0+6.1+0.7 

"^^■"^-7.6-13.0-0.5 

Q f;+3.3+3.4+10.8 

or; n+4.4+3.2+0.1 
Zd.U_44_g ]^_0.3 
PI q+0. 0+0.0+0.0 
"-'•"-'-0.0-0.0-0.0 
9 r,+0.0+0.7+22.8 
^•"-'-2.4-1.5-30.0 
/l-l O+17.5+13.0+24.3 
"^■•0-24.5-13.0-19.2 
90 Q+2.4+5.3+1.8 
^^■^-1.9-8.2-1.5 

q q + 0.4 + 0. 6 + 0.9 

'-••'-•-1.0-0.4-0.4 
n 7+0.2+0.2+0.6 
"^•' -0.2-0.0-0.3 

n no+0.01+0.00+0.55 

^■U^-0.02-0.01-0.52 


1.5 
-2.1 

-8.2 
-8.2 

13.3 



TABLE V: CP averaged branching ratios, direct CP asymmetries and the transverse polariza- 
tions fractions for Be I^*(s\^ decays, together with results from RCQM. 



g7^(10-^) 



7^r(%) 



channels 



Class 



This work 



RCQM 



This work 



RCQM 



This work 



Be 
Be 
Be 
Be 
Be 
Be 
Be 
Be 
Be 
Be 



D*+u 

D*+(j) 

D*+K*^ 
D*+uj 

Dl+o 



T 
C,A 
A,P 
A,P 
C,A 

P 

C,P 
A,P 
C,P 
A.P 



c;c: q+8.6+11.9+1.5 
00. n_ii.i 



5.0-11.1-1.4 
1.0+0.5+0.1 
0.8-0.6-0.1 
+59+5+11 
40-4-9 
+57+1+11 
42-1-9 
9 4+0.4+0.9+0.2 
^•^-0.6-0.7-0.1 

U.UU4_0-0-0.001 
+0.08+0.03+0.02 
-0.08-0.03-0.03 
^ 2+1.3+0.4+0.3 



3.^ 
I6I: 
172^ 



0.72 



-1.0-0.3-0.2 

+0.03+0.04+0.07 
0.01-0.05-0.04 
.-,+39, 3-0. 5+10. r, 
27.8-0.5-7.5 



0.26 



59.7 
13.0 
37.7 
30.6 



2.9 



38.8 



-24.1 



+3.0+4.2+0.4 



-3.4-2.7-0.4 
2Q 2+0.0+2.6+5.4 



^-1.5-5.8-7.6 



-14.9 



+1.1+3.1+0.3 
-0.8-1.7-0.1 



0.4: 



+0.0+0.0+0.0 
-0.0-0.1-0.0 
7 0+1.0+2.6+5.8 

" ' •°-0.0-3.4-5.0 



-1.3+7.6+1.4 
-1.1-4.5-0.9 
f. 9+0.1+1.3+0.0 
"^•^-0.3-1.8-0.1 

-5.3+6.8+7.9 

-4.6-6.8-4.7 
-0.1+0,1-0.0 
-0.1-0.1-0.0 



-29.3] 



-21. 3^ 
0.3^ 



3.8 
-3.0 
-6.2 
-0.8 



-3.0 
13.3 

-0.8 



-.p. /1+2. 5+2.0+0.3 

■■■ "3-^-1. 7-1.4-0.1 
54 3+1-8+4-0+0.5 



'-0. 9-2.4-0. 4 
+1.5+2.3+1.3 
-1.1-1.8-0.7 
+0.6+0.1+0.9 
0.7-0.1-0.4 
Q+1. 2+9.6+0. 7 

11.4 



52.6 
57.4 



-0.5-6.5-0.7 
22.3+0.0+5.3 
5.4-0.0-6.9 



1 1 9+0.5+2.1+0.2 
^-'-•'^-0.3 

CO Q + 2.1 

UO.O_2.3-4.4-0.4 



-1.4- 
-3.9- 



0.1 

0.8 



49.5 
67.E 



f8.8+2. 1+4.4 
-10.9-1.2-2.8 
+2.1+0,1-1.4 
-3.1-0.2-1.5 



for Be D+TT^. This makes Bn(Be D*+7r°) larger than Bn{Be D+n^). Simi- 
larly, we have Bn{Be ^ D*+p°) > Bn{Be D+p^). However, for Be D^*^+ri{ri'), 
while the dd part contributes to the annihilation diagrams, the constructive or destruc- 
tive interference situation are just the reverse, and BTZ{Be D*~^rj{rj')) are smaller than 
Bn{Be^ D+ri{ri')). 

For another kind of 6 — > s processes, the decays Be — >■ -Ds*^"^(7r°, a;) have a small 
branching ratio at O{10~^) due to the absent annihilation diagram contributions, and the 
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emission diagram contributions suppressed by CKM matrix elements Since the 

contribution of penguin operator is comparable to the one of tree operator, the interference 
between the two contributions is large. As a result, a big CP asymmetry is predicted in 
these decays. The branching ratio is even smaller ~ 10^^° and no CP violation for 
Be —7- D^*^'^(f) decays, since there are only penguin diagrams contributions. All these 
and other rare decays are also important, since they are quite sensitive to new physics 
contributions. 

For the ^*s)^ decays the branching ratios and the transverse polarization frac- 

tions TZt are given as 

where the helicity amplitudes Ai have the following relationships with A^'^''^ 

Ao = A"-, A± = A''± A^. (39) 

Our numerical results of the CP averaged branching ratios, direct CP asymmetries and 
the transverse polarization fractions for Be — t- D*_,^V decays are shown in Table IVl The 
transverse polarization contributions are usually suppressed by the factor ry or r£,* com- 
paring with the longitudinal polarization contributions, thus we do have relatively small 
transverse polarization factions for the tree-dominant decay {TZt{Bc — t- D*^p^) = 16.4%) 
and the pure penguin type decay {1Zt{Bc D*^(j)) = 11.5%). For the pure emission type 
decay Be — > D*'^u, the transverse polarization faction is large because the non-factorizable 
emission diagram induced by operate Oq can enhance the transverse polarization sizably. 
The fact that the non-factorizable contribution can give large transverse polarization 



contribution is also observed in the B^ — t- p^p^, uu decays 37|. For other decays, the 
annihilation type contributions dominate the branching ratios due to the large Wilson 
coefficients. Therefore, the transverse polarizations take a larger ratio in the branching 
ratios, which can reach 50% ~ 70%. This is similar to the case of i? — ?■ (pK* and various 



B pK* decays [38 



From Table |Vl one can also see that our branching ratios for B^ — )■ 
D*+K*^,D*^K*^,D*~^K*^,D*+(f) decays, are about 2 to 5 times larger than those in 
RCQM model, due to the sizable contributions of transverse polarization amplitudes. An- 
other point should be addressed that the annihilation contributions with a strong phase 
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have remarkable effects on tfie direct CP asymmetries in tfiese decays. As a result, our 
predictions are somewhat larger than those from RCQM. 

IV. CONCLUSION 

In this paper, we investigate the two body non-leptonic decays of Be meson with the 
final states involving one meson in the pQCD approach based on kx factorization. It 
is found that the non-factorizable emission and annihilation type diagrams are possible 
to give a large contribution, especially for those color suppressed modes and annihilation 
diagram dominant modes. All the branching ratios and CP asymmetry parameters are 
calculated and the ratios of the transverse polarization contributions in the Be — > D*s)V 
decays are estimated. Because of the different weak phase and strong phase from tree 
diagrams, penguin diagrams and annihilation diagrams, we predict a possible large direct 
CP-violation in some channels. We also find that the transverse polarization contributions 
in some channels, which mainly come from the non-factorizable emission diagrams or 
annihilation type diagrams, are large. 

Generally, our predictions for the branching ratios in the tree-dominant B^. decays are 
in good agreements with that of RCQM model. But we have much larger branching ratios 
in the color-suppressed, annihilation diagram dominant B^ decays, due to the included 
non-factorizable diagrams and annihilation type diagrams contributions. 
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Appendix A: factorization formulas for Be — > D*V 



We mark L, N and T to denote the contributions from longitudinal polarization, normal 
polarization and transverse polarization, respectively. 

J^e'"'' = 2^7rCffBfvM^ dx^J^ hb2dhdb2(pDl^{x2,b2)x 

{[x2 - 2rb + rDirb - 2x2)]as(ta)he{ae, Pa,bi,b2)St{x2) cxp[-Sab{ta)] 

+ [rl,{xi - l)]as{tb)he{ae, A, 62, bi)St{xi) exp[-Sabitb)]}, (Al) 

J-,"'^ = 2^7:CffBfvM%rv j\x2 b^b2db^db2ct>Di^p2M)>< 

{[rb-2 + rD{x2 + l- Arb)]as{ta)he{ae, /3a, bi, b2)St{x2) exp[-5'a6(ta) 
+rD[2xi - l\as{tb)he{ae, Pb, 62, bi)St{xi) exp[-Sab{tb)]}: (A2) 

•T^e"'^ = 2^7rCffBfvMl,rv 1^ dx^ h,h2dh^dh2(l>D;Jx2,b2) x 

{[n - 2 + roil - X2)\as{ta)K{,ae, /3a, bi, b2)St{x2) ex.p[-Sab{ta)] 
-rDas{tb)he{ae, /3b, 62, bi)St{xi) exp[-Sab{tb)]}, (A3) 

j^LR,L _ jrLL,L j^LR,N _ jrLL,N jrLR,T _ jrLL,T (A4) 

The factorizable emission topology contribution J-'^^'^[i = L,N,T) vanishes due to the 
conservation of charge parity 

8 /"^ 

M^^'^ = --TrCffsM^ dx2dx3 6263^^2^^^, (2:2, fe2)0y(a;3){[l - a:i 
3 Jo Jo 

-X3 - roixi + X2- l)]as{tc)he{/3c, a^, 63, 62) exp[-5'c<i(tc)] + [-1 + 

+X2 - X3 - roixi + X2- l)\as{td)K{^d, Oie, b3, &2) e^[-Scd{td)]}, (A5) 



8 

Me^'^ = -nCffsM^rv / dx2dx3 / 62MMMD; , (2^2, &2){[(2;i + X3 - 1) 
o Jo Jo 

(/)V{X3) + 2rD{x3 - X2)(t)v{x3)\0is{tc)he{P 63, 62) exp[-Scd{tc)] 
-[(-2r£,(l - 2xi -X2 + X3) - xi + X3)(l)\,{x3) + 2(r£,(l - X2 - X3) 
-2xi + 2x3)(f)v{x3)]as{td)he{/3d, cVe, b3, 62) ex.p[-Scd{td)]}, (A6) 
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8 /"^ 1"°° 

M^^''^ ^ -TiCffsM^^rv dx2dx3 h2hdh2dh?,^D' Xx2,h2) x 
-3 Jo Jo ^ ' 

{[{xi + X3- 1)0^(X3) - 2rD{2xi +X2 + XS- 2)0^(0:3)] 

Ols{tc)he{Pc: ^e, h: h) exp[-5'cd(ic)] 

-[(xs - a;i)(/.^(a;3) + 2{rD{2xi + X2 - X3 - 1) - 2xi + 2x3)(l)v{x3)] 
as{td)he{Pd: Q!e, ^3, ^2) exp[-5'c<i(t<i)]}, (A7) 

{[{Xi + X3-I+ rD{x2 - X3))(t)y{x3) + (Xi + ^3 - 1 - 

r£)(2xi + X2 + X3 - 2))(t)y{x3)\as{tc)he{l3c, ^e, ^3, ^2) exp[-5'cd(ic)] 
-[{xi - X3 - rc(l - X2 - a:;3))0y(a^3) - {xi - 2^3 + ^^dCI - 2xi 

-X2 + a;3))0y(a;3)]as(td)/ie(^d, Oe, 63, ^2) expf-^cdltd)]}, (AS) 

Q /•! /-oo 

^LR,T ^ ^LR,N ^_ (j^j^^A^^^ ^^^^^ / fe^MM^D^ (2:2, M 

63, 62) exp[-5'cd(tc)] 
+as{td)he{Pd: "e, &3, ^2) exp[-5'c<i(t<i)]}, (A9) 

8 /'^ f°° 

•^r'^ = -^^CffeM^ J dx2dx3 I h2hdh2dh(t>Di^^{x2M)<t>v{xz){[2-2x^ 

-X2 - X3 + roixi + X2- l)]as{tc)he{f3c, ae, h, 62) exp[-Scd{tc)] - 

[x3 - xi - roixi + X2- l)]as{td)he{Pd, ae, 63, 62) exp[-5'cd(td)]}, (AlO) 

O /"l poo 

M!''''' = --TrCffBM'^rv J dX2dXs J b2hdb2dh(f>D^Jx2,b2) X 

{[(xi + X3 - 1 - 2rD{2xi +X2 + X3- 2))0^(a;3) - (xi + 2:3 - l)0y(a;3)] 
asiQheiPc, ae, &3, ^2) exp[-5'cd(tc)] 

- X3){(j)y{xs) - (l)v{x3))as{td)he{Pd: ae, &3, ^2) ex.p[-Scd{td)]}, (All) 

Me^'^ = -TrCffsMsrv / dx2dx3 / 62&3C?&2(i&30Df , (a;2, ^2) x 
-3 Jo io 

{[(a^i + X3 - 1 - 2ri3(2xi + X2 + X3 - 2))(/)^(x3) - (xi + X3 - l)(t>v(x3)] 

as{tc)hei(3c, ae, &3, ^2) exp[-Scditc)] 

+ {xi - X3){(j)y{x3) - 4>v{x3))as{td)he{l3d, ae, 63, ^2) exp[-5'cd(td)]}, (A12) 
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/I poo 
dX2dX3 J b2hdb2db3(l)D*^^{x2,b2){[-X3(t)v{x3) 

+rcrv{(f>v{x3) - (f)y{x3))]as{te)he{aa, /3e, &2, ^3) exp[-Sef{te)]St{x3) + 

[X2(f)v{x3) + 2rvrDix2 - l)0y(x3)] 

as{tf)he{aa,/3f, 63, 62) exp[-Sef{tf)]St{x2)}, (A13) 

pi POO 

jrLL,N ^ -sCpfBTrM^rD J dx2dxs J b2bsdb2db3<PD'^^^{x2,b2){[rv{xs - l^^vi^s) 
-rc<f)v{x3) + rv{x3 + l)(l)v]as{te)he{aa, /3e, 62, ^3) exp[-Sef{te)]St{x3) 

-rv[{x2 + iWvi^s) - {X2 - l)rv{x3)] 

as{tf)he{aa, /?/, 63, ^2) exp[-Sef{tf)]St{x2)}, (A14) 



*' a 



/»1 /*oo 

= SCpfB-r^MBrD / dx2dx3 / b2b3db2db3(j)D*^^{x2,b2) x 
{[rvixs + 1)0^(X3) - re0^(x3) + rvix3 - l)(l>vi^3)] 



SP^L 



as{te)he{aa, /3e, ^2, ^s) exp[-S'e/(te)]5'f (^s) 

+ry[(-X2 - l)rv{x3) + {X2 - iWvi^s)] 

as{tf)he{aa,Pf, 63, ^2) exp[-,Se/(t/)]-St(a;2)}, (A15) 

■pLR,L _ jpLL,L jpLR,N _ jpLL,N ■pLR,T _ JpLL,T (AI6) 
pi poo 

IGCf/bTtMI / dX2dX3 I 62&3<^&2<^&30D; ,(3^2, &2){['^c0V'(a^3) 

Jo Jo ^ ' 

+rvX3{(f)v{x3) - (f)v{x3))]as{te)he{aa, Pe, &2, ^3) exp[-5'e/(te)] - 

[rDX20y(x3) - 2rv(l)v{x3)]ois{tf)he{aa, /3f, 63, 62) exp[-Sef{tf)]}, (A17) 

/I /-oo 
dx2dx3 J b2b3db2db3(t)D*^^{x2M) ^ 

{rD[4>v{^^) - '2,rcrv(t)v{x3)\as{te)he{aa, /3e, &2, ^3) exp[-S'e/(te)] 

+rv{(t>v{x3) + 0^(a;3))Qi.(i/)/ie(aa, /S/, &3, ^2) exp[-5e/(i/)]}, (A18) 

/I /-oo 
dX2rfX3 y b2b3db2db3(f>DlJx2,b2) X 

{rD[(pvi^3) - 2rcrv(f>v{x3)](^s{te)he{aa, /3e, 62, 63) exp[-S'e/(te)] 
+ry(0y(2:3) + 0y(x3))a.(t/)/ie(aa, b3, 62) exp[-5e/(t/)]}, (A19) 



*' a 



■SP,N 
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o 

M^^'^ = ^CfJbt^M^ I dx2dxs I bib2dbidb2(l)D*ix2,b2){[{xi-X2-rc) 

Jo Jo ^ ' 

4>v{x3) - rDrv[{x2 - X3)4>y{x3) - {2xi -X2- 2:3)0^(^3)]] 

as{tg)he{Pg, 61, 62) exp[- S gh{t g)] - [(1 -n-Xi- X3)(j)v{x3) 

-rnrvlixs - X2)(pv{x3) + (2x1 + X2 + X3 - 2)0^(0:3)]] 

as{th)he{Ph: C^a: ^1, ^2) exp[-Sgh{th)]}: (A20) 

16 f°° 
M^^'^ = -—CpfBT^MirDrv J dx2dx3 J hib2dh^dh2(t)Di^{x2M)(t>v{x^) x 

{rcO:s(tg)he{/3g, aa, bi, 62) exp[-Sgh(tg)] 

-rbas{th)he{/3h: oia, bi, 62) exp[-Sgh{th)]}: (A21) 

16 f°° 
M^^''^ ^ -—CpfBT^MlrDrv j dx2dx3 J bib2dbidb2(pDi^^{.X2M)4>vi.^i) ^ 

{rcas{tg)he{/3g, aa, bi, 62) exp[-Sgh{tg)] 

-nasithjheiPh: Ola: bi, 62) exp[-Sgh{th)]} , (A22) 



M^^^^ ^ -CpfBTrM^ dX2dX3 bib2dbidb2(pDiAx2,b2){[rD{xi-X2 + rc) 
o Jo Jo 

(pvixs) + rv{-xi +X3- rc)(0y(x3) + 0y(^3))] 

a,,{tg)he{/3g, aa, bi, 62) exp[-Sgh{tg)] - [-roixi + X2 - n - l)0y(a;3) 
+rv{xi + X3-rb- l)((/>y(a;3) + 4>v{x3))] 

Ols{th)he{Ph: Ola: bi, ^2) exp[-Sgh{th)]}, (A23) 

Q pi poo 

^LR,T ^ ^LR,N ^ _Cpfj,^M% J dX2dX3 J b,b2dbidb2cj>D*^^{x2, b2) X 

{[rv{xi -X3 + rc){(f)v{x3) + (Pvi^s)) - roixi -X2 + rc)0y(x3)] 

asitg)he{Pg, aa, 61, 62) exp[-Sghitg)] + 

[rv{xi +X3-rb- l){(f)v{x3) + 4>v{^3)) + ^^^(1 + n - Xi - X2)4'v{x3)] 
as{th)he{Ph, o^a, bi, 62) exp[-Sgh{th)]}, (A24) 
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8 /"°° 

M^^'^ = -CpfBT^M^ / dX2dX3 / bib2dbidb2(l)D* {X2,b2){[{xi - Xs - Tc) 

-J Jo Jo 

4>V{X3) - rorviix^ - X2)4>v{x3) - {2xi - X2 - X3)(py{x3)]] 

as{tg)he{/3g, aa, 61, 62) eicp[-Sgh{tg)] - [(1 -n-xi- X2)(f)v{x3) 

-rDrv[ix2 - X3)(f)y{xs) + {2xi + X2 + X3 - 2)0^(0:3)]] 
as{th)he{Ph: aa: ^1, ^2) exp[-Sgh{th)]}: (A25) 
_^sp,N^j^LL,N^ M^/'^ = -M^^'^. (A26) 

Appendix B: scales and related functions in hard kernel 

We show here the functions he, coming from the Fourier transform of hard kernel. 
he{a,P,bi,b2) = hi{a,bi) X h2{P,bi,b2), 

hiiaM) = It'i'^M:, """"^ 
^ ' ^ \ KQ{i^f^bi), a < 

hJR 6, h.) = [ ^(^1 - b2)Io{^b2)Ko{^b,) + {b, ^ 62), /3 > .3^. 

02) I ^^^^ _ i,^)j^^^i,^)K,{i^b,) + {b,^b2), /3 < ^""'^ 

where Jo is the Bessel function and Kq, Iq are modified Bessel function with Ko{ix) — 

^{—No{x) +iJo{x)). The hard scale t is chosen as the maximum of the virtuality of the 

internal momentum transition in the hard amplitudes, including l/bi{i — 1, 2, 3): 



where 



ta = max(v'|Q;e|, \/\Pa\, 1/^1, 1/^2), = max (^^, ^,1/61,1/62), 

tc = max(-v/|a^, a/I^, 1/62, 1/&3), = niax(A/|a^, a/}^, I/62, 1/&3), 

te = max(y|a^, \/|^, I/62, 1/&3), tf = max(v1aj, \/Wl^ 1/^2, 1/&3), 

= max( ^0^,-^^^,1/61,1/62), i/, = max(^4«J>^A^>l/^l>l/^2), (B2) 

Qie = {1-Xi- X2){xi - r^)M|, aa = -a;2a;3(l - rl,)Ml, 

Pa = [r^^ - ^2(1 - rl)]Ml, (3, = -(1 - x^){x^ - rl)Ml 

/3c = -(1 - - X2)[l - xi - X3(l - rl)]M],, 

Pd = {l-xi- X2)[xi -X3- rl{l - X3)]Ml, 

(3e = [rl - X3 - (1 - xs)rl]Ml = -^2(1 - r|,)]M|, 

= [^c - (^1 - ^-3(1 - rl)){x, - X2)]Ml, 

= [rl - (1 - xi - X3 + X3r^)(l -xi- X2)]Ml. (B3) 
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The Sudakov factors used in the text are defined by 



.Mfl , I Mb 



Sab{t) = 61) + 5(^0:2,62) + ^ / — 7g(/^)+2 / — 7g(/U), 

•J Jl/bi Jl/b2 



Mb Mb Mb Mb 

Scd{t) = s(--^Xi,62) ^X2,62) + S(--^X3,63) + -X3),63) 



4/ ^0,M + 2/ %M, 

Sef{t) = 5(^X2,62) + 5(^X3,63) + -X3), 63) 

+2 r ^7,(/i)+2 r ^7.(/i), 

Mb .Mb , . ^ .Mb Mb ^ 

-^xi, 61) + ^(^a;2, 62) + ^(^^3, &2) + 



+1/ ^7.(/i)+4/ ^7.(/i), (B4) 

where the functions s{Q, b) are defined in Appendix A of ^. 7g = —a^jT^ is the anomalous 
dimension of the quark. 



Appendix C: Light-Cone Distribution Amplitudes 

Here, we specify the hght-cone distribution amphtudes (LCDAs) for pseudoscalar and 
vector mesons. The expressions of twist-2 LCDAs are [3] 







-x)[l + 








-x)[l + 


4yC',/'{t)+alyCl/\t)], 




fT 


-x)[l + 


aiyC',^\t) + aiyCl^\t)], 



(CI) 

and those of twist-3 ones are 

= ^[l + (30r/3-^p|)C2^/'(t)-(r/3a;3 + ^P?.(l + 6ar))Cl/^(t)], 
Mx) = + 6(5^/3 - ^^3C^3 - ^pI - lpla^)il - lOx + lOx^)], 

where t = 2x — l, fy and are the decay constants of the vector meson with longitudinal 
and transverse polarization, respectively. For all pseudoscalar mesons, we choose rj^, = 
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0.015 and = —3 



26|]. The mass ratio p^^^x) = m^{K)liT\}^^ and p,,^^^^ = 2mq(^s)/mqq(ss), 



and the Gegenbauer polynomials C^(t) read 



cT{t) 



1 



3/2/ 



cTit) 
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'1 - Mr + 2ir 
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